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The locally observed cosmic ray spectrum has several puzzling features, such as the excess of 
positrons and antiprotons above 20 GeV and the discrepancy in the slopes of the spectra of 
cosmic ray protons and heavier nuclei in the TeV-PeV energy range. We show that these features 
are consistently explained by a nearby source which was active 2Myr ago and has injected 
(2 — 3) X 10®° erg in cosmic rays. The transient nature of the source and its overall energy budget 
point to the supernova origin of this local cosmic ray source. The age of the supernova suggests that 
the local cosmic ray injection was produced by the same supernova that has deposited °°Fe isotopes 
in the deep ocean crust. 


Introduction. Cosmic rays (CR) with energies at least 
up to 10^® eV are thought to be a by-product of the final 
stages of stellar evolution n The two main possi¬ 
bilities for the acceleration sites of CRs are individual 
supernovae (gamma-rw bursts, supernova remnants and 
pulsar wind nebulae) [l|, Q and superbubbles 0 hosting 
large numbers of supernovae (SN) and their progenitors, 
high-mass stars. 

The direct identification of CR sources which would 
allow the discrimination between these two possibilities 
is difficult, because the turbulent Galactic magnetic field 
(GMF) randomises the CR trajectories and leads to an 
almost isotropic CR intensity. Moreover, locally detected 
CRs are accumulated from a large number of sources 
which were active over the time scale Tesc 10 — 30Myr 
on which CRs (of the energy A ^ 10 GeV) escape from 
the Galaxy. The superposition of the signals from a large 
number of sources erases possible signatures of individual 
sources. 

The local CR flux might still have some ’’memory” of 
the individual sources composing it, because of the dis¬ 
crete and stochastic nature of the sources (be it SNe or 
superbubbles). The subset of near and recent CR sources 
could produce small features in the CR spectrum or cre¬ 
ate anisotropies The identification of such features 

could potentially provide a possibility for the identifica¬ 
tion of the CR sources and for the measurement of their 
characteristics. 

In what follows we show that the known differences in 
the slopes between CR protons and nuclei 9l-[lll . puz¬ 
zling features in the spectra of positrons [ia - [l^ and an¬ 
tiprotons mill could be self-consistently explained by 
a single nearby, recent CR source. We are able to de¬ 
duce the characteristics of the source from the details of 
the spectra of these CR flux components. In particu¬ 
lar, the hard spectra of antiprotons and positrons above 
^ 20 GeV and the soft spectrum of CR protons (com¬ 
pared to the spectra of heavy nuclei) can be explained by 
a source which has injected ^ 10®° erg in CRs in a tran¬ 
sient event which occurred ^ 2 million years ago. The 
source is located at a distance of (several) hundred par¬ 


secs along the local GMF direction. The transient nature 
of the event, its overall energy budget and the spectral 
characteristics of the injected GRs are consistent with a 
single SN and inconsistent with a superbubble as source. 

Contribution of a local source to the proton spectrum. 
Gosmic rays injected by a single source T years ago fill 
a region of the size dy _l ~ {D\\^±TY^‘^ in the inter¬ 
stellar medium (ISM). Here, and D± are the en¬ 
ergy dependent components of the diffusion tensor par¬ 
allel and perpendicular to the local GMF direction |1§ - 
[ 2 ^. If the total injected energy £tot is high enough, 
the source could produce a significant increase in the 
overall GR flux detectable by observers situated inside 
the region filled with CR. In the particular case of 
locally detected TeV CRs, the source contribution to 
the flux, F oc £tot/{d\\d\) oc could be com¬ 

parable to the locally observed CR flux if 
10®° erg/(l Myr®/^) - 10®^ erg/(10 Myr)®/^. Thus, 
both a SN which occurred a million years ago and has 
injected ~ 10®° erg in CRs and a super bubble which has 
injected ^ 10®^ erg over the last lOMyr can produce dis¬ 
tortions in the local CR spectrum. 

Cosmic rays spread faster along the direction of the 
GMF, ^ D±. The transient enhancement of the CR 
flux caused by a local source could be particularly strong 
if the source and the observer lie close to the same mag¬ 
netic field line. We model such a flux enhancement nu¬ 
merically using the code developed and tested in Ref. [ 23 . 
The code follows the trajectories of individual CR par¬ 
ticles through the GMF model of Jansson-Farrar 
starting from the moment of instantaneous injection in 
a single point by a transient GR source. The turbulent 
part of the field is chosen to follow isotropic Kolmogorov 
turbulence with the maximal length of the fluctuations 
Amax = 25 pc and the strength normalised to reproduce 
the observed B/C ratio, as discussed in [^J]. The 
calculation of the trajectories of individual CRs in the 
GMF allows us to include a detailed model for the regu¬ 
lar and the turbulent component of the GMF. We record 
the path length of CRs spent in a 50 pc sphere around 
the Earth, which can be converted to the local CR flux 
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at given time interval. 

We are interested in the case of a relatively young, 
T < fewMyr, and nearby source, dsource few x 100 pc, 
and CRs energies in the range 100 GeV-100 TeV. The 
spread of the CRs of such energy on Myr time scale is 
strongly anisotropic. A strong enhancement of the CR 
flux occurs if the source and the observer are connected 
by a magnetic field line. In this case, the contribution 
of a single source can dominate the observed total CR 
intensity at the Earth. 

Figure [T] shows an example of such a situation cal¬ 
culated for a source at the distance 300 pc which has 
injected CRs with spectrum dN/dE oc , 7p,inj = 

2 .2, and total injection energy Stot = 2.5 X 10®° erg. The 
source is placed at a CMF line passing within 50 pc from 
the Solar system. For A > 10 TeV, we calculate the CR 
trajectories up to 30 Myr, i.e. sufficiently long to observe 
the exponential cutoff in the flux due to CR escape. At 
any given energy, we find that the observed flux F at 
Earth as function of time rises, then drops as a power- 
law F{t) = up to the (energy-dependent) 

escape time and finally is exponentially suppressed as 
F{t) = exp(-t/Tesc). In the energy range 

1-10 TeV, we are only able to calculate trajectories up 
to 300 kyr. We extrapolate them to later times using the 
power-law with the slope a{E) derived from direct simu¬ 
lations in the energy range lOTeV-1 PeV. Note that the 
fluctuations visible especially at large times are due to 
the relatively small number of CR trajectories used. 

From Fig. [1] one can see that CRs with energies above 
100 TeV reach the Earth already 5 kyr after the injec¬ 
tion. If the source is able to accelerate CRs to energies 
above 10 PeV, their flux is suppressed already after 5 kyr, 
because of the fast escape from the Calactic disk. The 
escape induced flux suppression progresses towards lower 
energies with the increase of the source age. Below the 
high-energy cut-off, the slope of the spectrum softens and 
reaches the observed value 7 p ^ 2.7 — 2.8 after 2 Myr. 

The observed slopes of the spectra of the heavy nu¬ 
clei component of the CR flux, — 2.5, are systemat¬ 
ically harder than the slope of the proton spectrum in 
the TeV-PeV range [§, [T^. This harder slope of the nu¬ 
clear component of the CR flux consistently explains the 
shapes of the knees in the spectra of individual groups of 
nuclei within the escape model [1^ [2^ . The same slope 
of the average spectrum of O.l-lOTeV protons/nuclei in 
the Calaxy is deduced from a combination of gamma-ray 
and IceCube neutrino data [23 - 1^ . 

Assuming that the average Calactic CR proton flux 
at the Earth also has the slope jp ~ 2.5 and that 
it dominates the observed CR flux in the knee energy 
range E ^ 1-10 PeV (shown as “average” in Fig. [U see 
Ref. [13) one finds that the local source and the average 
Calactic contributions to the overall CR proton fluxes 
are comparable in the energy range 3-30 TeV. The uncer¬ 
tainty of the average flux is given by the uncertainty of 
measurements around the knee, cf. with the width of the 
blue band at PeV. The contribution of the local source 
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FIG. 1: Proton flux of the local source at different times. The 
average Galactic proton flux is shown as a thin orange line, the 
measured mectra of protons (light blue) from PAMELA [3, 
GREAM d, KASGADE and KASGADE-Grande [U as a 
band including experimental uncertainties. The sum of the 
average flux and the 2 Myr old source is shown by the blue 
thick line. 


with its softer spectrum explains the discrepancy between 
the slopes of the proton and heavy nuclei components of 
the TeV-PeV CR spectrum. In general, the local source 
gives also a contribution to the spectra of heavier nu¬ 
clei. If the elemental abundance of the local source CRs 
is identical to the overall measured CR abundance, the 
contribution of the source to the heavy nuclei spectra in 
the E > 1 TeV range is sub-dominant, because of the 
higher normalisation of the average Calactic component 
of the heavy nuclei fluxes. 

Positron excess from the local CR proton source. Our 
suggestion that the softer slope of the TeV-PeV pro¬ 
ton CR spectrum is caused by a local source can be 
tested via the identification of complementary signa¬ 
tures in the spectra of secondary particles—positrons and 
antiprotons—-produced in CR interactions in the ISM. 

The spectrum of CR positrons is known to have an 
’’excess” above 30 GeV. This excess refers to a devia¬ 
tion from reference models (as e.g. those of Ref. 113) 
which assume that positrons are solely secondary parti¬ 
cles produced during the propagation of time indepen¬ 
dent CR proton and nuclei flux through the ISM. The 
presence of this excess is usually considered as an indica¬ 
tion for the existence of a source of positrons in the local 
Galaxy. Source candidates under discussion are nearby 
pulsars [13, young (~ 10"^ yr) supernova remnants 
and dark matter annihilations or decays [30l| . 

A characteristic feature of secondary production in 
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FIG. 2: Spectum of positrons from the local source for the 
age 2 Myr (thin blue curve), 1 Myr (light grey) and 4 Myr 
(darker grey), compared to the measured spectra of positrons 
[ij-Tl^. The dashed orange line shows an estimate of the 
average Galactic positron flux [40l| . 


hadronic interactions is that the slope of the energy spec¬ 
trum of secondaries is very close to the one of the parent 
protons, since scaling violations are small except close 
to mass thresholds. Thus, the presence of a local source 
of CR protons should reveal itself through an associated 
component in the CR positron spectrum. The average 
energy fraction transferred to positrons produced in a CR 
interaction is {ze+) ~ 3% for 7 = 2.2, so that positrons 
with energies in the 30-300 GeV range are produced by 
protons with energies > 1-10 TeV. The flux of this lo¬ 
cal positron component is a function of time. It grows 
as the fraction of protons interacting with the ISM in¬ 
creases, but keeps approximately the shape of the parent 
proton distribution. Our calculation of the CR trajecto¬ 
ries emitted from the local source presented in the previ¬ 
ous section allows us to determine the average grammage 
X traversed by the CRs since the moment of injection. 
For energies E < 10^"* eV, the grammage is nearly en¬ 
ergy independent, X ce 0.3g/cm^, for a source of the age 
T = 2 Myr. 

The produced positrons diffuse and spread over larger 
and larger distances, softening thereby their energy spec¬ 
trum. The process of diffusion and the resulting softening 
of the spectral slope is identical for positrons and protons, 
if the age of the local source is small enough that energy 
losses of the positrons can be neglected. Thus, not only 
the injection, but also the propagated spectra of protons 
and positrons from the local source have nearly identical 
slopes at any moment of time. In particular, this implies 
that at present 7 e+ ~ 7 p ~ 2.7 — 2.8. 

Figure [ 2 ] shows the measured positron spectrum in the 
30-300 GeV energy range |12l - [l^ together with our cal¬ 
culation of the positron flux from the local source. For 
the calculation of the hadronic production cross sections, 
we have been employing QGSJET-II-04 [3l| in the mod¬ 
ified version presented in |32| . Since we have fixed both 
the contribution of the local source to the proton flux and 
the grammage in Sec. 1, the normalisation of the shown 


positron flux is a prediction. Both the normalisation and 
the slope 7 e+ — % agree well with the experimental data. 

An additional suppression of the positron 

flux may occur due to synchrotron and in¬ 

verse Compton energy losses. The synchrotron 
and inverse Compton cooling rate is = 

0.5[17/(0.5 eV/cm^)][Ae+/(300 GeV)]Myr“^, where 

U is the combined energy density of radiation and 
the magnetic field, U = 0.5 [i3/(4 ^G)]^ eV/cm^. The 
synchrotron and inverse Compton cooling softens the 
positron spectrum. Contrary to the primary CRs, 
which are injected instantaneously from a point source, 
positrons are continuously produced at a constant rate. 
This leads to the formation of a cooling break from 7 e+ 
to 7 e+ -I- 1 in the positron spectrum. The break energy 
decreases with time as shown schematically in Fig. [51 
The non-observation of such a softening limits the age 
of the local source to T 4 Myr. 

Note that the association of the observed positron ex¬ 
cess with a local source also implies a lower limit on the 
source age, T > 2 Myr. Otherwise, if the source would 
be much more recent, CRs would not have enough time 
to produce the observed excess positron flux. 

Overall, the hypothesis of the local source contribu¬ 
tion to the CR proton spectrum passes the positron self- 
consistency check and provides an explanation to the ob¬ 
served excess in the positron spectrum. 

Antiproton flux from the local CR source. Interactions of 
CR protons with the ISM also produce antiprotons. The 
relative size of the antiproton and positron fluxes is in the 
regime of negligible positron energy losses completely de¬ 
termined by ratio of the corresponding Z-factors, or ap¬ 
proximately by the ratio of the spectrally averaged en ergy 
fraction {zi) transferred to antiprotons and positrons [3^. 
We use again the modified version of QGSJET-II-04 
for the calculation of the spectrum of secondary antipro¬ 
tons, keeping all parameters fixed to the ones used for 
positrons. The resulting spectrum of secondary antipro¬ 
tons from the local source is shown in Fig. |5|by the blue 
line. The light-blue shaded range shows the ±50% uncer¬ 
tainty band attributed to the uncertainty of the antipro¬ 
ton production cross-section in this energy range (32i - l^ . 

From this figure one can see that the flux of the local 
source antiprotons constitutes a significant fraction of the 
overall GR antiproton flux. The local source contribution 
to the antiproton to proton ratio is at the level ^ 10 “"^, 
which is comparable to the measured ratio [T1,[T3. Note 
that the very preliminary AMS-02 antiproton measure¬ 
ments reported in 11 Til are consistent with the published 
PAMELA results [l6l | up to 180 GeV. The apparent inde¬ 
pendence of the pjp ratio with energy is at tension with 
the naive expectation that in the steady state model of 
secondary antiproton production by CR interactions in 
the ISM the antiproton spectrum should be softer than 
the proton spectrum, so that the ratio should have de¬ 
creasing trend with the increase of energy (35j | . The pres¬ 
ence of the single source contribution to the antiproton 
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FIG. 3: Spectrum of antiprotons from the local source (dashed 
blue line) compared to the measurements from [l^ . [l7| (thick 
light blue line). The orange box shows an estimate of the aver¬ 
age Galactic antiproton flux from [40)|: the horizontally orange 
hatched range shows an alternative calculation from [^. The 
blue hatched ranges show the uncertainty of the model calcu¬ 
lation (inclined hatching) and for the sum of the local source 
and average antiproton fluxes (vertical hatching). 


flux removes this tension. The energy dependence of the 
p/p ratio for the single source and for the average Galac¬ 
tic CR flux are different. Antiprotons originating from 
the source are injected by the parent protons which have 
spent the same time in the ISM, independently of their 
energy. For this component, the p/p ratio is rather slowly 
rising for energies well above the antiproton production 
threshold. The resulting independence of the p/p ratio 
with energy is therefore a falsifiable prediction of our 
model. Overall, the presence of the local source com¬ 
ponent in the antiproton flux is consistent and possibly 
even favoured by the p/p ratio data. 

Supernova nature of the local CR source. The com¬ 
bination of the positron, antiproton and proton signa¬ 
tures of a local CR source provides a possibility to con¬ 
strain its parameters. The source should have the age of 
T ^ 2Myr. An older source would be inconsistent with 
the absence of radiative cooling in the positron spectrum. 
A younger source would fail to produce sufficient amount 
of antimatter. A younger source would provide also a 
harder feature in the proton spectrum in the TeV range, 
while for an older source the source contribution would 
be too soft to be noticeable in the TeV band. 

The overall energy injected in CRs should be at the 
level of £tot 10 ®° erg, for a wide range of source dis¬ 
tances. The energy density of high-energy particles in¬ 
side a region in the ISM filled with CRs injected by the 
source is nearly uniform. The region spans a ~ 100 pc 
wide tube of kpc-scale length the along the local GMF 
direction. The only condition for the detectability of the 
CR flux from the local source is that the Solar system is 
situated inside this region. 

The source has to be transient, in the sense that it 
could not remain active all through the last million years. 
Otherwise the source would strongly contribute to the 


CR flux from more recent injection times. For exam¬ 
ple, the contribution from the last 0.1 Myr period would 
have the spectrum shown by the orange curve in Fig. 
[T] multiplied by O.lMyr/lMyr = 0.1. This contribu¬ 
tion would be comparable to the 2 Myr old contribution. 
However, the energy spectrum of this younger contribu¬ 
tion is harder, with the slope 7 ~ 2.5. The presence of 
such a contribution would erase the effect of softening 
and the spectrum of TeV-PeV protons would have the 
same slope as the spectrum of heavy nuclei. 

The transient nature of the source and the overall in¬ 
jected energy rule out the possibility that the local source 
is a superbubble blown by massive star formation. The 
typical lifetime of superbubbles is in the 10 ^- 10 ® yr range 
determined by the lifetime of massive stars. A young su¬ 
perbubble formed a million years ago would be still active 
today. 

The only plausible model of the transient local CR 
source is that of a SN. The average rate of SN explosions 
in the Milky Way disk volume Vdisk — 100-300 kpc® is 
Ti-SN — (1“3) X 10“^yr“^, or one SN per (0.3-3)x 10^yr 
per kpc®. This means that one could reasonably expect 
that about one SN has exploded within the last million 
years within a 100 pc wide, kpc long filament directed 
along the GMF line going through the Solar system. 

Discussion. Our analysis has shown that several fea¬ 
tures in the CR spectrum which appear puzzling within 
the standard Galactic CR injection/propagation models 
And a natural explanation by the presence of a local CR 
source. In particular, we have shown that a ^ 2 Myr old 
source which injected ~ 10®° erg in CRs with energies 
up to at least 30 TeV can consistently explain the differ¬ 
ence in the slopes of the proton and heavy nuclei spectra 
in the TeV-PeV energy range, give an additional con¬ 
tribution to the antiproton spectrum in agreement with 
recent AMS-02 data and predict the correct amplitude 
and slope for the positron spectrum in the 30-300 GeV 
energy range. The local source also gives rise to an excess 
anisotropy of the CR spectrum in the 1-100 TeV energy 
range, which for first time explains the anisotropy data 
in this energy range [4l|. Note that the combination of 
the anisotropy and the positron data constrain the source 
age and disfavor suggestions of younger sources in denser 
environments as in Ref. (i^ . 

The scenario of a 2 Myr old local source contributing to 
the 1-100 TeV CR spectrum could be tested with future 
AMS-02 and ISS-CREAM measurements of heavy nuclei 
CR fluxes. Such measurements could detect the imprint 
of the local source on the primary-to-secondary CR nuclei 
ratios. In particular, the energy independent grammage 
(~ 0.3 g/cm^) traversed by Carbon nuclei should result 
in the presence of an energy-independent component of 
the B/C ratio in the 1-100 TeV energy band. The 2 Myr 
time scale is also comparable to the decay time of ^°Be 
nuclei. This may leave an imprint in the ^°Be/°Be ratio 
in the same energy range. 

The presence of a nearby SN explosion was previously 
noticed in a completely different type of data on abun- 
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dance of isotopes on Earth [Sdl - fs^ . These data suggest 
that an episode of deposition of ®°Fe isotopes in the mil¬ 
lion years old deep ocean crust was produced by the pas¬ 
sage of an expanding shell of a 2 Myr old supernova rem¬ 
nant through the Solar System. The consistency of the 


SN age and distance estimate from the deep ocean sedi¬ 
ment data with those found from CR data suggests it is 
one and the same supernova which is responsible for the 
CR injection and the isotope deposition on the Earth. 
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Appendix A: Supplementary material 
1. Trajectory approach 

In the standard approach to Galactic CR propagation, 
the iV-particle phase space distribution of CRs is approx¬ 
imated as a macroscopic fluid. Solving then the diffusion 
equation including CR interactions permits to model self- 
consistently primary and secondary fluxes. However, the 
diffusion approach has several disadvantages. First of all, 
the diffusion approximation itself breaks down in several 
interesting regimes [H-Q- Moreover, the diffusion tensor 
is an external input in this approach. It is only loosely 
connected to fundamental properties of the Galactic mag¬ 
netic field (GMF). In this work, we apply therefore the 
trajectory approach where one calculates the path of in¬ 
dividual CRs solving the equations of motion of particles 
propagating in the GMF. We use the code described and 
tested in y, while we use for the regular and the tur¬ 
bulent component of the GMF the parameters described 
in Q. As we showed previously in the resulting 

CR fluxes reproduce over a wide range of energies all 
available experimental data for individual groups of CR 
nuclei. 

Our numerical results for the propagation of CRs emit¬ 
ted by a bursting source can be summarized as fol¬ 
lows. In an initial period, CRs propagate essentially 
one-dimensional, filling a filament-like structure around 
the magnetic field line going through the source. Af¬ 
ter this rather short initial period, CRs start to diffuse 
also perpendicular to the direction of the filament. In 
a realistic GMF model, the two perpendicular compo¬ 
nents of the diffusion tensor differ, D±^ 2 , and 

thus CRs diffuse mainly in the || and the _L, I direc¬ 
tion. In particular, we found for the case at hand in the 
Jansson-Farrar model @ the values D±^i, D±^ 2 } 
{10^°, 10^®, 10^^} cm^/s for CRs with energy lOOTeV. 
This behavior is illustrated in Fig. 0] which shows the 
scaled CR density n at the position of the Earth as func¬ 
tion of time for three different energies. Additionally, we 
used the relation D{E)t oc valid for Kolmogorov 

turbulence to rescale the Eq = 10^'^ eV curve by t —>■ f = 
to A = 3 X IQi^eV and to A = 1 x lO^^eV. 
The approximate agreement of the two rescaled curves 
with the true result supports a Kolmogorov-like scaling 
of the density. Moreover, it is visible that the density 
scales approximately as n oc 1/t for intermediates times, 
lOkyr <t< lOMyr. Such a scaling is characteristic for 
two-dimensional diffusion and reflects the observed hier¬ 
archy in the eigenvalues of the diffusion tensor. 

2. Secondary fluxes 

The trajectories of charged CRs depend on their rigid¬ 
ity TZ = E/{Ze). The sign difference of the charge of 
protons and antiprotons becomes important only on dis¬ 
tances of a few coherence lengths to the source, while on 



Time [kyr] 


FIG. 4: The number density n{E,t) of CRs at the position of 
the Earth as function of time for three different energies; also 
shown are the densities rescaled using i = (Eo/E)^^^ (grey 
lines) and the 1/t behavior expected for d = 2 diffusion. 

larger scales the diffusion of protons and antiprotons of 
the same energy proceeds in the same way. In this paper 
we do not perform the full (computationally expansive) 
calculation of trajectories of the secondary particles. In¬ 
stead, we use a scaling relation to describe the spread of 
the secondaries, taking into account the shift of energy 
of the secondaries compared to the primary particles. 

The density np{E', r',t') of antiprotons which are pro¬ 
duced at the time t by protons with the density np{E, r, t) 
is given by 

np{E',r',t') = J d'''x Gp{r',t'-r,t)np{E,r,t), (Al) 

where Gp{r 'is the Green function describing the 
diffusion of antiprotons. Our numerical results show that 
the CR propagation from a bursting source can be mod¬ 
eled by two-dimensional diffusion in the time and energy 
range of interest for CR secondary production. Thus the 
proton density at the production time t is given by 

exp [-rVriiff] , (A2) 

’’diff 

where d = 2, r denotes the distance in the two- 
dimensional diffusion plane, Q{E) the source spectrum 
and Tdiff the effective diffusion distance, r\g = ADt. The 
corresponding Green function G has the same functional 
form, 

Gp{r ' ^ exp [-(r - r'f/r'ls] (A3) 
^diff 

with r'^ig = AD'{t' — t). For simplicity, we assume a 
one-to-one relation between the energy of protons and 
antiprotons, setting E' = a?E ~ O.IE, where a ~ 0.46. 
This results in the diffusion coefficient D' = aD for Kol¬ 
mogorov turbulence. 
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The integral in Eq. mil is of the standard Gaus¬ 
sian type. In the limit of small r the exponent is 
exp [—~ 1 and the ratio n-^jn^ is 


R 


V — 


% 

rip 


1 

X -|- a(l — x) ’ 


(A4) 


where we introduced also x = t'/t. In particular, we see 
that the ratio does not depend on energy and satisfies as 
required Rp{a = 1) = Rp{x = 1) = 1. Integrating over 
the production time t, we find {Rp{x)) = ln(a)/(a — 1) ~ 
1.43 for a ~ 0.46. Thus, using for CR antiprotons the pri¬ 
mary proton trajectories we underestimate the antipro¬ 
ton flux by ~ 50%. Since this enhancement is energy 
independent, we can use the following scheme for the 
calculation of the antiproton flux. 

In a first step, we calculate the proton flux at the 
position of the Earth neglecting interactions. In the 
trajectory approach, we obtain the local CR flux at a 
given time interval recording the path length of CRs 
spent in a sphere around this point. We record also 
the grammage X these CRs crossed on their way to 
the Earth. Then we calculate the average grammage 
(X) = I dt p{xi{t)) summing up the gas den¬ 

sity along the trajectories of the N individual CRs. We 
employ n{z) = p/rup = no exp(—( 2 ;/zi/ 2 )^) as model 
for the gas distribution in the Galactic disk, with z as 
the distance to the Galactic plane, no = 0.3/cm^ at 
the Solar position and Zi /2 — 0.21 kpc Q- We set also 
n = 10“^g/cm^ as minimum gas density up to the edge 
of the Milky Way at j^l = 10 kpc. Next we calculate the 
secondary antiprotons produced by CRs interacting with 
gas using the modified version of QGSJET-II-04 intro¬ 
duced in 0. The antiproton flux is then given by 


dpi^Rp) — Rp^nucXjlTlp Zp(^Ep^cRj Ipi^E^^ (A5) 


where the Z-factor is given by the inclusive spectra of 
antiprotons dap{E, zp)/dzp, Zp = Ep/E, 


Zpi^Ep^od^ — / dz > 
Jo 


dap{Ep/z,z) 

dz 


(A 6 ) 


with a ~ 2.6 as the average slope of the proton spectrum. 
The factor Rp ~ 1.5 accounts for the underestimate of 
the antiproton flux, while we include the effect of heav¬ 
ier nuclei in the interstellar medium using a nuclear en¬ 
hancement factor enuc 1-2. Note that it is a virtue of 
our model that the slope of the injection spectrum is the 
same for all CR nuclei. As a result, the use of a nuclear 
enhancement factor is justified. 


The same approach could be applied also to positrons, 
with a change a ~ 0.31 to take into account their lower 
average energy. The above simple picture has to be modi¬ 
fied at the highest energies (beyond the current measure¬ 
ments energy range) where continuous energy losses of 
positrons become important. The first estimate of the 
corrections due to the continuous energy losses could be 
found from the solution of the diffusion equation in the 
Syrovatskii approach in d = 2 M, 


n{E, r) 


Q{E)b{Ei) 


exp[-r7r^ig] , 


(A7) 


where r^jg = 4 J*dtD{E), b{E) = —dE/dt = /3E‘^ the 
ener gy l oss rate and Ei denotes the energy at injec¬ 
tion [^. The effective diffusion coefficient for this so¬ 
lution is 13 ^ Z3|| as determined from our simulations. 
For a fixed value of (3, this determines the electron and 
positron spectra formed at time t. The total electron and 
positron fluxes are then obtained summing over the for¬ 
mation time. This approximation has to be verified by 
simulations of trajectories of secondary positrons taking 
into account continuous energy losses. Such simulations 
require however a considerable increase of computing re¬ 
sources and are therefore postponed to a future work. 
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